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BRIEF 


Water samples and geological materials are analyzed for as many 
as 28 elements by using 0.8 to 4,0 MeV charged particles to excite 
X-rays and ^-rays which are monitored simultaneously. 



abstract 


Charged particles ranging in energy from 0.8 to 4.0 MeV are used to 
induce resonant nuclear reactions. Coulomb excitation (2f” ra y s )> an< 3 X-ray 
emission in both thick and thin targets. Quantitative analysis is possible 
for elements from Li to Pb in complex environmental samples, although the 
matrix can severely reduce the sensitivity. It is necessary to use a 
comparator technique for the $ -rays while for X-rays an internal standard 
can be used, A USGS standard rock is analyzed for a total of 28 elements. 
Water samples can be analyzed either by nebulizing the sample doped with 
Cs or Y onto a thin formvar film or by extracting the sample (with or without 
an internal standard) onto ion exchange resin which is pressed into a pellet. 



INTRODUCTION 


The possibility of using protons accelerated to energies of several 
MeV as an excitation source for X-rays or y-rays has attracted much 
attention. The bombardment of a target by heavy charged particles 
frequently results in the ejection of inner shell electrons and the 
subsequent emission of characteristic X-rays* By using lithium drifted 
silicon, SI (Li), detectors for energy dispersive analysis, the emitted 
X-rays can be used for qualitative and quantitative analysis 
of any element heavier than Na. in a relatively simple sample as little 
as 10~ 12 gm of an element can be detected (l). In addition to the 
interaction with electrons, accelerated ions take part in a variety of 
nuclear scattering and reaction processes. In the MeV energy range, 
the cross sections for these reactions exhibit sharp resonances associated 
with the formation of metastable compound systems between the projectile 
and the target nucleus with a lifetime long compared with typical 
periods of internal nuclear motion. Protons are especially good for 
inducing these resonant reactions with light nuclei. By controlling 
the energy of the bombarding protons it is possible to selectively analyze 
for one element in a complex matrix because of the narrowness of the reson- 
ances., By increasing the energy of the protons above the resonance 
energy to compensate for the energy loss with penetration depth, this 
technique has been used as a nuclear microprobe to show depth distribution 
of the element of interest (2). The y-rays emitted in these nuclear 
reactions are characteristic of the compound or product nucleus. 
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Another type of nuclear excitation is due to the interaction of the 

population 

Coulomb fields of the projectile and target nuclei, resulting ir^of excited 
nuclear states with subsequent decay by y-ray emission* The dominant 
mode of excitation is through the electric quadrupole (E2) interaction 
(3), thus only nuclei with appropriate low lying nuclear states will 
be excited in this manner* 

Each of the three types of interactions has advantages for analytical 

applications. The atomic excitation of X-rays is sensitive for elements 

heavier than N% but complex matrices can limit sensitivities* Resonant- 

nuclear excitation is selective and sensitive for light elements even 

in complex samples. The Coulomb excitation gives specific y-rays for 

nuclei of each isotope but is much less sensitive. A complete 

analysis of a target is possible by simultaneously measuring both X-rays 

and y-rays. By using the energetically well defined beam from a Van de 

Graaff accelerator the light elements are measured one at a time in the 

order of increasing resonance energy* The X-rays are measured at one 

proton energy. The greatest problem with this type of analysis is In the 

integrated 

conversion of /\ i peak areas for X-rays and y-rays to accurate precise 
elemental concentrations. This paper will describe procedures used for 
elemental analysis of complex samples such as biological and geological 
materials B by^ a combination of the proton excitation processes described 


above , 
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theoretical 

In the bombardment of a target with a charged particle beam, the 

probability of a nuclear or atomic reaction is expressed as a cross- 

section which is a function of the kinetic energy of the beam. When the 

electromagnetic radiation (X-rays and y-rays) emitted in these reactions 

is monitored, the yield (Y . » events/unit charge) or count rate in a 

particular detector will depend upon the total cross section (c?(E) * 

cm 2 /atom) for formation of the excited state at the bombarding energy 

E, the product (k *)Tfk.-) of the attenuation factors (k-) 

for the emitted radiation 

/or the fraction of excited nuclear or atomic states which result in 
emission of the radiation being observed, the density of target nuclei 
(N m particles/cm 3 ) in the target material, the thickness (X = cm) of 
the target, the inverse of the ionization charge state of the particles 
in the beam (l = beam particles/unit charge), the efficiency (c) of the 
detector for the radiation energy being observed, and the fractional 
solid angle (_TL) subtended by the detector, 

dY * NI-flekcr(E)dx (l) 

Charged particles readily lose energy in penetrating the target. 
This energy loss can be expressed as the areal stopping power S(E) which 
is the rate of energy loss per unit of path length divided by the total 
target density (/^ “ * gm/cm 3 where is the partial density of 

element i in the target). 

c dE 


S(E) 


( 2 ) 
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Thus , if the average energy of the particles in the beam impinging on the 
target is E q , the energy of the particle in the target is between E q and 
the lower average energy (E^) of the emergent beam which has penetrated 
the thickness (X) of the target. For a thick enough target the beam is 
stopped and E^ = 0. Since the cross section is a function of the energy, 
equation 1 must be integrated over the energy range from E q to E^. 
Integrating equation 1 and substituting for dx a dE/(dE/dx) gives: 


Y = NI _TL.sk 


E. rr(E)dE 


\E q (dE/dx) 


(3) 


N 

Since N “y^i — (where N - Avogadro's number and A. = gram atomic 
i 

weight of element i), using the areal stopping power gives the yield for 
element i in the target; 



(4) 


The fraction of the partial densities due to element i = c^) is 
the concentration of i in the target. In the case of X-ray excitation, 
the energy resolution obtainable with a Si(Li) detector is not sufficient 
to measure differences between the isotopes of element i and the elemental 
concentration can be extracted by measuring the characteristic X-ray 
yield ^Y?j according to equation 4. 

In the case of nuclear reactions it is necessary to consider the 
fractional abundance ^a_ ^ of the jth isotope of element i. 
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The yield for the particular reaction becomes: 


Y. . = ifie.k 

xj A , J x 

I J I 


^.(EjdE) 


c .. a. , 
L 


A convenient simplification may be introduced where thin targets 

/V E f \ 

(j-g «1 ) are employed. If J E q — e £ | * AE is sufficiently small such 

that one may consider cr (E ) - constant over the domain of &E > then, 


o(E)dE ^ 

"stll 


(e)t 


where T «* £> X 
_ X 

energy at 


areal target thickness, and E ® the average beam 


'uantitative Analysis 


Absolute or Comparison Methods 


Absolute determination of elemental concentrations requires a good 
knowledge of the areal target thickness or, in the case of thick targets, 
the areal stopping power function S(e), For complex matrices, this 
information is not often easy or even practically possible to obtain. 

If the target matrix consists of one or two very dominant elements 
which are known, such that the function S(E) may be considered essen- 
tially to be determined by them, then analysis by absolute methods becomes 
practical. However, this is not necessarily the general case. 
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Comparison methods consist of three basic types, (a)those where 
the target sample is doped with a known amount of a known element, 

(b) those where the concentration of a certain element in a sample is 
determined by some other method and subsequent measurements of other 
elements are made relative to it, and (c) n comparator" methods ( 4 ) 
where elemental or isotopic ratios may be determined by reference to 
a second prepared standard containing the elements or isotopes sought 
in known relative abundance. 


Methods (a) and (b) are identical in principle. Using the subscript 
R to refer to the reference element or isotope, one may manipulate 
eq. 5 to obtains 


a 

a 


c 

1 1 


R C R 



\Vr 

V k 

K 1 1 


a R (E)dE 
E q S(E) 

7 (E)dE 
1 . 

E S(E) 

O 


(7) 


If, in addition, a thin target is used, one may employ eq. 6 to arrive 
at the following: 


a c 
_J A. 
a„Ci 
R R 



k R e R A i?R^ 

riTA-fm 

1 1 R 1 


( 8 ) 


Little additional error is added if is substituted for E, provided 
the targets are sufficiently thin. However, for thick samples one still 
requires a knowledge of S(E), except in certain instances. If, 
for example, we write a = where B is a constant with respect to 
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bombarding energy, then the ratio of integrals in eq. 7 may be replaced 
by the quantity B -1 ; this is the case, at least approximately, for 
certain instances of Coulomb excitation and X-ray production. 

A formula for the stopping power (S(E)) proposed by Bethe and later 
modified by Bloch which is valid for high energies and a single element 
is given by: 


S i (E) 


4TTe 4 N 

J Q 

M c ? A, 
o i 





(9) 


Here, e * electron charge, M q * electron mass, c = speed of light in 

vacuo, N q = Avogadro's number, A^ * atomic weight of stopping material, 

Z. * atomic number of stopping material, Z, = atomic number of projectile, 

V i /eT7 

0 * — s .0463 / 7 M:. where E. is the bombarding particle energy (MeV) 

C y 1 1. 

and is the particle mass (amu), and I is a quantity related to the ioniza- 
tion potential of the stopping medium and is given by: 

2 / 

1^9.1 Z t [1 + 1.9 2 ' 3 ] eV (10) 

In principle, the stopping power for a complex matrix can be derived from 
this, using 

S(E) =^Tc.S.(E) ( 11 ) 

where the c, are the elemental concentrations (weight percent). Experimen- 
tally determined values for SjE) have been compiled in a review by North- 
cliff and Schilling (5). 
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In the "comparator” methods, it is necessary for the total stopping 

powers of the unknown and standard target samples to be identical within 

certain error limits. The samples and standards need not be exactly 

amounts 

the same but similar in the relative/\of major elements as shown by Ricci ( 4 ) 

for biological materials and by Sipple and Glover (6) for sedimentary 

rocks. Hence, letting primed quantities indicate the standard sample, 

and recognizing that the same radiation is being observed by the same 

detector in the same configuration for both the standard and unknown 

and assuming equality of the stopping powers (S^S*)^ then equation 7 

becomes: v 

ac Y 

£ T 7 r "7 T ( 12 ) 

To be completely valid, the beam of charged particles should be completely 
stopped in both samples; otherwise, the difference in target thickness 
must be considered. 


experimental 

In this work a 5.5 MeV Van de Graaff accelerator is used to produce 

beams of protons, alphas, or neon ions at well defined energies. Beam 

currents ranging from IQ. nanoamps to about 2 microamps are 

focused and collimated on the targets. The y-rays from resonant nuclear 

reactions and Coulomb excitation are , measured using an Ortec Ge(Li) 

detector. A Kevex Si (Li) detector is used to measure the X-rays from the 
excitation. 

atomic A Information from each detector :is fed through a separate ADC 
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into one half of the memory of a ND3300 (Nuclear Data) multichannel 

analyzer. The 7-ray spec trim in one half of the memory and the X-ray 

are • - 

spectrum in the other half/\read out on magnetic tape for computer 

analysis on a CDC 6400 at a later date. The experimental arrangement 

is shown in Fig* 1. During bombardment the target chamber as a whole 

current. 

is electrically insulated to allow integration of the beam/\ The Si(Li) 
detector could be included in the vacuum system or separated by a 
thin (25.4/t(m) Be window* The Si(Li) detector insures good efficiency 
for X-rays but is insensitive to the higher energy p otons resulting 
from nuclear reactions* In the case of some y-rays emitted in reactions 
the energies are so high that the Ge(Li) detector is not very efficient. 

In these cases greater sensitivity is possible by using a large volume 
Nal(Tl) detector since the y-ray resolution is not very important for 
the resonant reactions where usually only one type of nuclide is excited. 

The simple target holder allows up to 12 samples to be analyzed without 

breaking the vacuum (and also holds a quartz viewer for beam focussing). 

Both thick and thin targets are used for a large variety of samples 

including biological and geological materials. Thick targets are pre- 

3000 psi 

pared from dry powders by pressing at A in a die in a hydraulic press. 

When there is not enough adhesion, samples are mixed with graphite 

(Crescent Dental Mfg.) 

powder with a Wig-L-Bug^before pressing. The resulting pellets ? which 
were 1.1 cm in diameter are attached to 1.9 cm A1 squares with a small 
amount of silver conductive paint (manufactured by C.G. Electronics for 
printed circuit boards) or epoxy and these targets are mounted onto the 
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sample holder* Thin slices of fine grained rock are" also used* Several 

types of thin targets are prepared by evaporating solutions onto. 

various backings* Drops of solutions placed on Ta or A1 disks give non- 

uniform residues upon evaporation. In order to make the residue uniform, 

a nebulizer is used to spray the solution onto the backing material* 

The method described by Jolley and White (7 ) is used to reduce solid 

biological and geological samples by ultrasonif ication to the liquid 

form, allowing them to be sprayed onto the backing with the nebulizer 

Model No. 

(DeVtlbiss Nebulizer y^lSO) . In addition to the thick metal backings, 
thin p Id stic backings are mounted on an Al washer 1.9 cm in diameter 
with a center hole 1.1 cm in diameter* The samples are then sprayed 
onto these thin films of formvar, mylar, polystyrene, or carbon. Commer- 
cially available films are greater than 6.35/^m thick but can be made 
thinner in the laboratory. The advantage of these thin low Z backings is 
that the background they produce by bremsstrahlung or excitation is 
negligible* 

After a series of samples \ar.e placed in the target chamber the 

energy of the beam is set at the lowest energy of the resonant 

reaction to be measured and then focused onto a quartz viewer. In 

succession, each sample is, bombarded for a fixed integrated charge 

the 

without changing the beam and/\re suiting y-ray spectra are recorded. Then 
the beam energy is increased to the next resonance energy and the process 
repeated for each sample. At the highest resonance energy the X-ray 
spectra are recorded along with the Y^ray .spectra. 
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RESULTS AND DISCUSSION 

"X-ray Excitation" 

The main advantage of using heavy charged particles for X-ray 
excitation is that the background can be quite low while the cross 
section for X-ray production remains high* The background associated 
with the heavy ion Is bremsstrahlung caused by the secondary electrons. 
Classically the maximum energy imparted to an electron in a head on 
collision can be calculated as four times the product of the ratio of 
electron to ion mass and the energy of the ion. Thus a 4 MeV alpha 

particle , which can excite X-rays above 20 keV, produces electrons 
with less than 2 keV of energy, keeping the bremsstrahlung background 
confined to low energies. This can be compared to using electron 

excitation where the bremsstrahlung covers the entire energy region 

to 

being excited or/\excitation with monochromatic X-rays which because 
of Compton scattering has the highest background just below the energy 
of the exciting line, Cairns (s) sites an improvement in the peak/ 
background ratio of from 10 3 for 20 keV electrons to 10 5 for 2 MeV 
proton excitation with comparable X-ray yields. The lower energy of 
the background for heavy ion excitation can be easily suppressed by 
simple filters ( 8). However, the use of absorbers adds complications to 
accurate quantitative analysis. When thin targets are used, the back- 
ground is, of course, much smaller than for thick targets. 
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Figure 2 compares the X-ray spectrum for proton excitation on a 
thick graphite-containing target (a) and a thin film of formvar (b). 

In both cases the sample is the USGS standard rock BCR-1. In (a) 

0.1984 gm of the rock is mixed with 10$ graphite powder (National Spectro- 
scopic grade SP-1C) which acts as a binder and also provides electrical 
conduction, 

before being pressed into a pellet. In the thin sample (b) 0.002 gm 
of BCR-1 is nebulized onto a lOO/^g/cm 2 thick formvar film. Proton 
current on the sample .is imited to 25-30 nA$ and a total of 40 .micro- 
Coulombs of charge is collected^ since only the major elements are 
sought. Note the difference in th$ 3 to 5 keV region under the 
X-ray peaks for K, Ca, Ti and V, A 20 mg/cm 2 A1 absorber is placed 
between the thick sample (a) and the Si(Li) detector to decrease the 
total count rate by eliminating some of the low energy bremsstrahlung. 

Several types of thin targets are prepared from formvar (lOO^gm/cm 2 ) , 
mylar (640/{g/cm 2 ),,and polyvinylidene chloride (l200Xfgm/cm 2 ) and analyzed 
for background. In all of these materials the largest amount of back- 
ground is in the low energy (2-3 ke\^ bremsstrahlung region. The 
intensity of the background increases with increasing target thickness. 
Obviously the characteristic Cl X-rays in polyvinylidene chloride make • 
it uridesireable , The formvar is superior to the mylar because of 
its higher heat resistance allowing use of higher currents (up 

to 100 nA with no cooling). The formvar films prepared in this laboratory 
are thinner than the commercial mylar used which means that the 
background is nearly two orders of magnitude lower for the formvar. 

The .one advantage of mylar is that the position of the beam on the 
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target could be observed due to the fluorescence in the mylar. For 

the 

both the thin hydrocarbon backings andK thick graphite- containing pellets 

an additional source of low energy background is from the periodic 

discharge which occurs when a positive charge builds up on the nonconducting 

targets. This is a significant source of background in the thin target 

experiments but it can be largely eliminated by placing a bare filament 

in the target chamber to spray therm 3 ! electrons onto the target, thus 

preventing positive charge accumulation ( 9 )« This background source 

containing 

is eliminated from the thick graphite-fltargets by attaching them to the 

the 

Al backing with the conducting silver paint rather than^epoxy. It 
is not possible to assign absolute values to the sensitivity obtainable 
by this technique as it depends upon the experimental arrangements, , e.g. 
the detector solid angle and other characteristics. However, the 
relative sensitivities can be quite well described by the production 
cross section (the ionization cross section times the fluorescence yield). 

Experimental measurements (10, ll) of cross sections suggest that 
the binary encounter approximation (BEA) predicts the ionization cross 
sections very closely. The fluorescence yields have been compiled 
by Bambynek et al. (12). Figure 3 shows the ionization cross sections 
for the three types of particles used in this study (H + , He + , Ne + ) as 
a function of the atomic number of the target. These measurements 
were made using the University of Virginia 5.5 MeV Van de Graaff 
accelerator (10). In general this indicates that protons 



-14- 


are the best particles for excitation because for a given particle 

energy the ;ionization cross sections are greater by at least an order of mag 

nitude^ However, in complex samples there may be a problem with too 

many X-rays excited for the detector to function well. In this case 

the heavier particles (like neon) are more selective in that elements 

within rather limited ranges of atomic numbers are efficiently excited. 

In addition Cairns (8) has pointed out that below a certain critical 
energy a bombarding heavy ion will not excite X-rays, Thus, by careful 
choice of heavy bombarding ion and energy, a particular element can be 
selectively excited. For a more complete characterization of the sample, 
however, the protons are more universal. In addition to the cross 
sections, the background differences and depth of penetration must 
also be considered. The bremsstrahlung background Is at lower energies 
and the depth of penetration is much lower for the heavier particles. 

An absolute measure of the sensitivity is determined by nebulizing 
and evaporating 10 ml of a dilute aqueous solution onto a formvar 
film and measuring the yield. Table 1 indicates these sensitivities 

obtained by extrapolating the peak area obtained down to a value of 

obtained 

twice the uncertainty in the background from the same integrated 
(1000^ C) 

charge on a blank formvar film. These calculations are based on a 
solid angle of 5,5 x 10” 4 steradians subtended by the Si(Li) detector. 
Another method of analyzing for trace elements in aqueous solutions is 
to extract onto an ion exchange resin which can be pressed into a pellet 
(using 10$ graphite). Table I also gives the sensitivities obtained by 



15- 


extracting 10 ml of the same solutions as used for the thin films 
onto 0,250 gms Bio-Rad Chelex 100 resin 0 The advantage of this tech- 
nique is that larger liquid volumes could be rapidly batch extracted 
onto the resin, thus improving the sensitivity. 

The greater sensitivity obtained by using the ion exchange resin 
also 

is due to problems with the nebulization technique. The 

spray covers a larger area than the resin pellet <s but the size of 

Therefore, 

the proton beam is constant- a smaller fraction of the evaporated target 
is sampled. In addition^ our experiments have indicated that there is 
also substantial loss of the material in solution due to deposition 

inside the nebulizer. The pellets give a uniform and reproducible 

-it he 

target 1.1 cm in diameter. In contrast the sample on thin film, 
limited by the A1 washer to the same size, is not uniformly distributed 
and hence is less reproducible. Thus in terms of 

both time required and greater sensitivity, the resin technique is best 
for routine water analysis. For more complex samples such as geological 
materials, analysis is more difficult and less sensitive. This is due 
to the high intensity of the K X-rays for the relatively light major 
elements. These obscure the low intensity L X-rays in the same region, 
cause pulse pile up in the electronic circuitry, and degrade the 


resolution. 
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Two methods of quantifying the analysis of aqueous solutions are 
investigated. In the comparator method a standard solution is used to 
prepare a target containing the elements of interest. The X-ray intensi- 
ties from the sample can be compared directly with those from the standard. 
This does require that conditions be kept uniform for the sample and 
standard. For instance, the targets must be uniform and reproducible 
and either the areal density or total weight of the sample must be known. 
The non-uniformity of the thin film targets prepared by nebulizing and 
evaporating makes precision impossible using the comparator technique. 

The only precise way to analyze the thin targets onto which aqueous 
solutions have been nebulized is to use an internal standard technique. 

In the case of the resin pellets the precision is very good. Individual 
pellets, each containing Cs and three other elements which were extracted 
from 1 ml of solution, are prepared. Another pellet is prepared by 
extracting from a mixture of the nine elements. This pellet is used as 
a standard to calculate the concentrations in the other resin pellets. 

Table II shows the results of using the comparator method for solutions 
extracted onto the resin. The accuracy is good (within expected uncer- 
tainties in counting statistics) except in cases where selective absorption 
and enhancement effects occurred in the standard mixture (not shown). 

The necessary corrections for this type of effect can be made using 
techniques such as those described by Rasberry and Heindrich (13), 
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sections , the 

By knowing the production cross /\ mass absorption coefficients, 

and the thickness of any absorbers between the sample and detector, 

absolute concentrations can be determined by comparing the intensity 

e lement 

of the X-rays from each element to an/\which has been added to the 

original solution at a known concentration, Gove et al. (14) suggests 

CsN 0 3 as the doping material since it is soluble and the ions do not 

form precipitates. Both the 4,29 keV Lc\and 30,8 keV K^X-rays can 

be used. In this work, and Gove*s (1.4), the cesium concentrations are 

varied (100-1000 ppm) to ensure that the results were independent of 

the dopant concentration. The standard solution used for the comparator 

method is : doped with Cs and the concentrations measured by the internal 

standard method are shown In Table II. In this case 5 ml of the doped 

solution is nebulized onto the formvar backing. Relatively poor 

precision and accuracy can result from the high background (due to 

the charge buildup on the target) under the Cs L X-ray, The background 

could be eliminated with thernal electrons (.9) or could be avoided 
Yittrium for Y 

by using .. as the internal standard. The X-ray energies^are Kc^ (14.93) 

and L^(l.92)« and it is not a common contaminant in water samples, is 

soluble as the nitrate, and does not interfere with any of the other 

elements which might be sought in environmental samples. The resin 

samples could also be analyzed by this internal standard technique. 

The results for the mixture of nine elements is also shown in Table II. 

It should be noted that the uncertainties are higher for the internal 

standard results than for the comparator results due to ; the uncertain- 

cross 

ties in production4 sect £°hS. Again corrections for selective enhancement 
and absorption in the thick samples are necessary. 
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’’Nuclear Resonant Reactions” 

The nuclear reactions used are either inelas.tic proton scattering 
with subsequent Jf"-ray emission (p,p r ^O or proton absorption followed 
by emission of an alpha particle and/or decay to the ground state 
( P» %). ' The cross sections for these reactions 

vary with energy, being small except for narrow energy bands or resonances 
where they are as much as several orders of magnitude higher. The 
energies and widths of the resonances may be found in the nuclear 1 
physics literature. One example is shown in Fig, 4 for three £7 A1 i 

i 

(p,^) 28 Si resonances, A thin Al film (~100 / ^g/cm a ) is vacuum 
deposited on a 0,25 mm Ta backing. The energy of the protons bombarding 
the target is increased in steps of 2 keV and the counts for 300 micro- 
Coulombs of integrated charge are recorded. Note the slight tailing 
of the peaks on the high energy side due to the finite thickness of 

the Al layer. In fact, since the rate of energy loss for the 

(dE/dx) 

proton/^ can be calculated, it is possible to use these resonant 

of 

reactions as a measuremeny|depth (2). For instance , deposition of another 
layer of material (other than Al) on the Al surface changes Figure 4 
by shifting the peak up on the energy scale to a degree which depends 
on the thickness of the covering layer. 

At the resonance energies the reactions listed in Table III are 
selective in the sense that only one reaction is prevalent, however * 
at higher energies (especially for thick samples) other reactions 



-19- 


also occur. The reactions used must be chosen to yield unique products 
(a unique $-ray spectrum) to avoid interferences between different 
nuclear reactions yielding the same final isotope. The unique ray 
spectra of the products and the enhanced cross sections at resonance 
energies combine to make this a very selective technique. Sensitivities 
are measured using thick target pellets formed from USGS standard 
rocks mixed with graphite. Thin targets are prepared by nebulizing 
solution onto thin formvar films 9 but the resulting non-uniform target 
does not allow the necessary precision. Vacuum deposition is used 
to put a uniform thin layer of the desired element on a Ta backing 
(as in Fig. 4). The minimum detectable limit is defined as a ^-ray 
peak (minus background) larger than two standard deviations of the back- 
ground in that region of the spectrum. The results are tabulated in 
Table III. 

In the case of the thin films the thickness is known so that 
sensitivities can be calculated in units of mg/cm 2 . If the area 
covered by the beam is known (in this case measured to be approxi- 
mately 0,07 cm 2 ) sensitivity can be calculated in terms of grams on 
the surface of > the target. When the proton beam is set at the resonance 
energy the thickness of material sampled is very small and depends 
upon the width of the resonance. Thus s in the case of the thick 
geological samples, the sensitivities (in ppm) may appear poorer because 
such a small portion of the sample is actually being analyzed. 



- 20 - 


For quantitative analysis the comparator method must be used. Since 
Van de Graaff accelerators can be held at a particular energy a series 
of samples should be analyzed for a single element before going to the 
proton energy required for the next element. The only critical factor, 
besides the energy, which must be kept constant is the angle of the 
sample relative to the beam and the ray detector. It should be noted 
that this nuclear microprobe technique can be used to measure spatial 
distribution with a well focussed beam (15) or to measure the average 
concentration using a defocussed beam as was used in this work. 


"Coulomb Excitation" 


The excitation of the nucleus by the interaction of its Coulomb 

field with that of the bombarding nucleus is a purely electromagnetic 

process. The cross section for Coulomb excitation depends on several 

factors including the energy, spin and parity of the excited nuclear 

state being populated relative to the ground state, the energy of the 

bombarding particle, and the strength of the Coulomb interaction. 

a Z Z 2 

The last is described by the Sommerfeld number F a (Z and 

— * nv p 

Z^_ are the nuclear charges on the projectile and target nuclei respectively 
and V is the velocity of the bombarding particle). The classical treat- 
ment of this phenomena (3) indicates that unless fe^‘4 use of this excitation 
for analysis would not be practical. This relationship means that Coulomb 
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excitation is more probable for targets with high atomic numbers and 
for the heavier alpha bombarding particles. This, however, does not 
take into account the fact that the nucleus must have a low lying 
excited state. As the energy of this excited state increases, the 
excitation probability decreases, but at different rates for protons 
and for alphas. The ratio of the cross sections using 4 MeV protons 
versus 4 MeV alphas is calculated as a function of the excitation 
energy for four different values of Z^. and plotted in Figure 5, Thus, 
the decision as to the best type of bombarding particle must be made 
with respect to the energy level being excited. It should be kept 
in mind that the proton has a greater range in a thick target and, 
therefore, samples a larger volume. However, even under the most 
favorable circumstances the cross sections for these reactions are low. 

In Table IV the sensitivities for a few of the most favorable 
isotopes are given. These are measured by analysis of thin films of 
the elements vacuum deposited on Ta backings. Again the vacuum deposition 
is used because the nebulized aqueous solution is not as uniform. In 
addition the Ta backing allowed greater beam currents and thus more 
rapid data collection. The limit of 0.1 counts^^ Coulomb for a 
Ge(Li) detector 10 cm from the target is set on the basis of the low 
background for these samples (except for Ta peaks). In more complex 
samples, such as geological materials (Fig. 6), the background is higher 
due to the increased Compton scattering of the -rays produced 
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by some of the resonant nuclear reactions listed in Table III. This 
increased background decreases the sensitivity by as much as two 
orders of magnitude. 

If the Coulomb excited ^-rays are used to analyze samples, a 

comparator technique is necessary. Although the sensitivities are not 

. and 

very good, this phenomenon is quite selective^ so could be of use for 
particular applications such as measuring one of the elements in Table 
IV in a matrix which does not yield Coulomb excited ^-rays. In this 
work it is found that some of these elements could be measured in 
samples at the same time as the X-rays and the y«rays from the nuclear 
resonant reactions. Thus, while perhaps not a valuable analytical 
technique by itself, when combined with the other two types of excitation 
it can provide useful information. 

Analysis of Geological Materials 

For solid materials, e’.g. rocks or lyophilized biological samples, 
the best sensitivity is obtained by forming a pellet using 10# graphite. 
These samples can withstand larger beam currents than the thin films 

i 

although the thick target makes it necessary to correct for the loss 
' bombarding 

in energy of the^|particle as it passes through the sample. The stopping 
power of a thick sample can be considered as a matrix effect since 
different volumes of material are sampled. If the internal standard 
method is used, then the only matrix effect which need be considered is 
the differential absorption of X-rays of different energies which can be 
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taken into account using programs such as those described by Brown et al. 
(l6)« In order to use the comparator method the mass stopping power (l/^p) 
(dE/dx) must be known. Since most geological materials are composed of 
a few major elements in a narrow range of atomic numbers, the stopping 
powers do not vary appreciable between different samples. Table V gives 
the values of ( 1/jO ) (dE/dx) for five different proton energies in a 
number of USGS rock standards (17) calculated using the data of Northcliff 
and Schilling (5). Although these igneous and metamorphic rocks cover 
a wide range of compositions, the stopping powers for a particular energy 
are very similar. For each energy the average and the deviation from 
the mean are given. This indicates that by using one of these standardized 
rock powders for comparison, almost any type of rock can be analyzed without 
a correction for the small differences in stopping powers. It should 
be pointed out that the volume of specimen actually sampled is small 
since the depth sampled depends upon the range of the bombarding particle. 
The proton beam used for analysis is on the order of 3 mm in diameter 
and the actual amount of rock pellet being sampled by a 4 MeV proton is 
calculated (5) to be on the order of 3,5 x 10“ 3 gms. Therefore, unless 
the sample is carefully powdered and mixed, sample inhomogenieties may 
lead to imprecise results. In the case of the nuclear resonant reactions 
an even smaller amount of specimen is effectively sampled, the actual 
depth depending 'upon - ' the width of the resonance. 
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The proton resonance and X-ray analysis of replicate samples of 
the same IISGS standard are reproducible within the uncertainties 

due to counting statistics. The most critical factor is found to be 
controlling the angle of the target relative to the beam and detector. 
The necessary precision should be built into the sample holder. 

In order to show the potential of a combination of the three 
techniques described above, several USGS standard rocks are analyzed, 

A granite (G~2) is used as the primary standard and the concentrations 
are calculated by comparison. The results of one of these rocks (a 
basalt, BCR-l) are given in Table VI along with the range of values 
reported by Flanagan (17), Some of the deviations from the best 
values reported for BCR-1 may be due to sample inhomogenie ties on the 
macro or micro scale and some (such as for Li) are thought to be partially 
due to problems in proper positioning of the sample. The data for the 
other rocks is as good p but not as complete, as the BCR-1 results. 

The analysis sequence for a set of samples such as these geological 

t 

standards is as follows: l)prepare pellets using ~10$ graphite as 

binder* 2):mount samples with silver conductive paint onto Al squares, 
3)align the samples on the holder and mount in sample chamber, 4)after 
evacuating the chamber set proton energy at 0,991 MeV(Bs resonance) and 
focus the beam on the quartz window, 5)i rradiate the standard (G-2) 
with 0,5^amps and count the ft -rays for l,000^_Coulombs or until adequate 
counts are obtained (as long as the collected charge is measured), 
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6) storey) -ray spectrum on magnetic tape and move the first sample into 
the beam, 7) repeat data collection for each sample at a particular 
energy and then change to energy of next resonant reaction and analyze 
the standard and each sample, 8)after finishing the P resonance at 
3.73 MeV increase proton energy to 4.0 MeV and collect both X-ray and 
^ -ray spectra for SO^Coulombs . An example of the X-ray spectra is 
that for BCR-1 shown in Figure 2a. The ^-ray spectrum for BCR-1 obtained 
at the same time is shown in Figure 6. 

The photopeaks in the ^-ray spectra taken at 4 MeV are due mostly 
to nuclear reactions. Since the target is thick at different depths 
within the sample the proton has lost enough energy so that it is at 
one of the resonance energies for each reaction. The intensities of 
these peaks are used to calculate the concentrations of the different 
elements and the results agree with those taken earlier for the analysis 
of the surface. This indicate^ that the samples are homogeneous. In 
addition to these large photopeaks there are several peaks due to 
Coulomb excitation. The peaks observed are noted in Table VI along 
with the calculated concentrations using G=2 as the standard. The lack 
of sensitivity is readily apparent. The concentrations obtained for 
the light elements using the nuclear resonant reactions are also given. 

In the case of the X-rays, corrections are tqa.de for absorption and 
secondary fluorescence by some technique such as that of Brown et al. (16). 
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This type of correction requires the analysis of more than one standard 

to initially calculate intensity ratios. Once these are calculated a 

single standard can be used as a monitor for each set of unknown samples. 

could be improved by 

■ The sensitivity^longer counting and irradiation times. 

Use of the energy dispersive Si(Li) detector for the analysis of 

complex samples suffers from the lack of resolution, A detector with 

that 

higher resolution thanAused for this work would eliminate some of the 

overlap of peaks. In order to reduce the resolution problem a deconvolu- 

program 

tion program was written for the CDC 6400. This^fits the background 
region of the spectrum with an exponential or power series polynomial 
and subtracts the fit from the spectrum. The spectrum with background 
subtracted is examined for photopeaks using the standard deviation 
above background and numerical first and second derivatives at each 
data point as criteria. Multiplets found from the maxima and minima of 
the second derivative are resolved with the option of fitting each 
multiplet to a function consisting of N gaussians. The program provides 
peak limits, centroid, and area, for each resolved peak. 

The analysis of these geological samples is an indication of the 
possibilities for analyzing complex samples, A total of 28 elements 
are measured simultaneously. With longer irradiations and a better 
Si(Li) detector more elements could be measured. The greatest 

sensitivities are, however, obtained by analyzing solutions nebulized 
and evaporated onto thin formvar films or extracted onto ion exchange 

resin which is formed into a pellet. As discussed above, the X-ray 
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intensities can be converted to concentrations using the internal standard 
or comparator method. The ^-rays from resonant reactions or Coulomb 
excitation can be measured at the same time as the X-rays, but quantitative 
analysis requires a standard containing the elements of interest for 
comparison. Preparation of thin uniform targets and the limitations in 
beam current are problems. This is a very sensitive technique for 
environmental studies or liquid samples. Although a large particle 
accelerator is necessary for this type of analysis, time on these machines 
is becoming increasingly available for analytical work. The calculations 
can be more straightforward than for neutron activation analysis if the 
samples are prepared as described above so that the rather complex 
equation 5 can be simplified. 
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TABLE I. 

Detection Limits for 4.0 MeV Charged 

a 

Particle ^Induced X-tays 

Element 

FORMVAR 
Protons 
(10~ 9 gm) 

Alphas 

(10 -9 gm) 

RESIN 
Protons 
(KT 9 gm) 

Cl 

K a 

12.2 


1.40 

Co 

Ka 

5.9 

7.6 

0.81 

Ni 

K a 

8.3 

11.2 

1.18 

Ge 

Ka 

15.2 

72.5 

1.33 

Se 

Ka 

46.8 

75.0 

1.67 

Rb 

Ka 

52.1 

246.0 

2.11 

Cd 

Ka 

165 

-- 

9.13 

Cs 

Ka 

-- 

-- 

79.5 

T1 

La 

42.3 

109 

3.33 


Limits defined for the minimum detected by nebulizing (for 
formvar) or extracting (for resin) 10 ml of solution and 
measuring X-rays for 1 ,000^Coulombs of integrated charge. 




TABLE XI. 

Determination of Elemental Concentrations 




Using 

4.0 MeV Proton- 

-Induced X-rays 



Si 

Cesium Internal Standard 

Comparator 

Method* 3 


Element 

Formvar 
Backing, ppm 

Ion Exchange 
Re s in, ppm 

Ion Exchange 1 ' 
Resin, ppm 

Actual, ppm 

Cl 

17,6 + 

2,5 

16.2 + 

2,4 

. — , 

16.3 

Co 

27.4 + 

3.6 

27.6 + 

3.6 

25.9 + 0.2 

25.3 

Ni 

31,1 + 

4.0 

26,8 + 

3.4 

25.3 + 0.3 

27.0 

Ge 

24.6 + 

3,0 

22.8 + 

2.8 

w fc* 

25.6 

Se 

26.9 + 

3.2 

26,8 + 

3,1 

- 23,1 + 0,2 

25.8 

Rb 

30.7 + 

3.8 

27,6 + 

3.2 

— 

24.7 

Cd 

29 ,9 + 

3,9 

22.6 + 

2.5 

24.9 + 1.8 

26.0 

Cs 



— 


232 + 17 

250 

Tl 

25,1 + 

2,0 

26.1 + 

2.1 


25.8 


Solutions doped with 250 ppm Cesium, 

b Mixture of all nine elements used as the standard. 



£ 

TABLE III. Detection Limits for Resonant Nuclear Reactions 


Reaction 

Proton 

Energy 

Gamma 

(MeV ) Energy (MeV) 

Thin target 
Limit ,/<gm 

Thick targ* 
Limit , ppm 

7 Lt (p, p'#) 7 Li 

1.030 

0,478 

0.10 

0,69 

9 Be(p,2) 10 B 

0.991 

7.50 

-- 


10 B(p,#) ll C 

1.532 

0.429 

0,84 

<2200 

11 B (p , p l ^0 1I 'B 

2.664 

2.144 

0.70 

<2000 

12 c(p,^ 13 n 

1.700 

3.51 

— 

«* m 

ls N(p,a3) 12 C 

0.896 

4.430 


-- 

19 F(p, a33 1S 0 

1.375 

6.130, 7.120 

1.3 

0.96 

23 Na(p, p'$0 23 Na 

1.458 

0.439 

0.024 

0.58 

24 Mg(p,p'<30 24 Mg 

2.930 

1,370 


28.7 

27 Al(p s p , 2T) 27 Al 

2.727 

1.013 

— 

23.3 

as Si(p ,p'20 2s Si 

3,100 

1.779 

510 

62.3 

31 p(p,pW s1 p 

3 .730 

1.266 

12.0 

1*1 

32 S(p,p*#) 32 S 

3.379 

2.237 

1.8 

3.5 

3S Cl(p,p'# 3S Cl 

2.721 

1.220, 1.763 

2.1 

6.2 

a Based on 1000^ Coulombs of integrated charge. 



L 

°Thin target (100/^g/cm 2 ) prepared by vacuum deposition of material onto 

0.25 mm thick Ta 

backing. 




c Thick targets (l 

mm) were 

prepared by pressing 

a mixture of rock powder 


and 10# graphite. 



TABLE IV ♦ Detection Limits for 4.0 MeV Charged Particle 
Induced Coulomb Excitation Gamma Rays 


Element 

Isotopic 

Abundance 

Gamma 

Energy* keV 

Minimum 

Weight 

Target, 

detectable 
for Thin 
As, 

Minimum detectable 
Concentration in 
Thick Target, mg/gm 

47 Tt 

7.3 

159 

.. , J 1 

proton 

0.48 


alpha 

0.15 

proton, 

0.2 

alpha 

0.7 

4Srji ^ 

73.9 

983 

0.44 


6.30 

0.2 

30.0 

ss Mn 

100 

126 

0.03 


0.01 

0.01 

0.04 

58 Fe 

91.7 

847 

0.32 


3.78 

0.1 

18.0 

57- 

Fe 

2.2 

136 

2.17 


0.68 

1.0 

3.2 

75 As 

100 

199 

0.26 


0.10 

0.002 

0.5 

77 Se 

7.6 

239 

0.77 


0.38 

0.3 

1.8 

so Se 

49 ,8 

666 

0.46 


4,13 

0.2 

19.7 

lsl Ta 

100 

136 

0.08 


0.04 

0.03 

0.2 

isi Ta 

100 

301 

0.48 


2.24 

0.01 

10.7 

197 Au 

100 

192 

2,52 


3.22 

1.2 

15.3 

197 Au 

100 

279 

0.51 


2.31 

0.2 

11.0 

a 

Based on 

0.1 cts/^Coulomb detected 

In Ge(Li) 

10 

cm from 

target. 




table 

V. Stopping Powers 

for USGS Rock Standards 3 


Proton 
Energy, MeV 

G-2 

PCC-1 

USGS STANDARD 
GSP-1 

BCR-1 

G-l 

1.0 

194.76 

194.95 

193.71 

189.49 

194.94 

2.0 

119.40 

119.41 

118.89 

116.75 

119.50 

3.0 

89.04 

89.05 

88.71 

87.29 

89.11 

4.0 

71.80 

71.81 

71.56 

70.50 

71.85 

5.0 

60.72 

60.73 

60.53 

59.68 

60.77 


Proton 
Energy, MeV 

DTS-1 

AGV-1 

W-l 

Average 

Standard 

Deviation 

1.0 

195.25 

192.66 

190.13 

193.24 

2.28 

2.0 

119.56 

118.35 

117.09 

118.62 

1.13 

3.0 

89.14 

88.35 

87.52 

88.53 

0.74 

4.0 

71.88 

71.29 

70.67 

71.42 

0.55 

5.0 

60,79 

60.31 

59.82 

60.42 

0.44 

Stopping powers 
(17). 

in units 

of keV*cm s /mg, 

based on 

major element 

compositioi 



TABLE VI. Measurement of Elemental Concentrations 
in Geological Standard BCR~1 & 


Element 

Reported f ppm 

Proton 

Resonance „nr>m 

X-ray 

Excitation „ppm 

Coulomb 

Excitation .ppm 

Li 

12.8 

12.4 + 0.8 



F 

485 

491 + 9 

— 


Na 

2.42* 

2.46 + 0.03$ 

— 

— 

Mg 

, 2.08* 

2.06 + 0*04 $ 


— 

A1 

7.20* 

7.12 + 0.2 .$ 

— 

— 

Si 

23.82* 

23 .8 + 0.3 $ 

— 

— 

P 

0.157* 

0.140 + 0.03 $ 

— 

warn 

S 

392 

384 + 10 


-- 

Cl 

50 

43.2 +2.5 

tm o* 

m» m* 

K 

... 1.41* 


1.10 + 0.52 $ 

~ 

Ca 

4.94* 

— 

4.37 + .0.30 $ 

— 

Ti 

1.27* 

— 

1*24 + 0.02 $ 

1.13 +0.28 $ 

V 

399 

— 

390 + 4 

334 + 72 

Cr 

17.6 


18*5 + 0*5 

— 

Mn 

1406 


1390 + 12 

1238 + 240 

Fe 

9.37* 

— 

9.14 + 0* 02 $ 

8.46 + 1.52$ 

Cu 

18.4 

““ 

17.7 + 0.6 

— 

Zn 

120 

— 

118 + 2 

106 + 28 

Ga 

20 


19*5 + 2*2 

— 

As 

0.70 


0*70 + 0*02 


Rb 

46.6 


45.2 + 1.2 

— 

Sr 

330 


332 + 0.4 


Y 

37,1 

Itnft *Xft 

38.2 +2.1 

— 

Zr 

190 


188 + 2,4 

— 

Nb 

13,5 


15*4 + 2*0 


A g 

0,036 


0.032 + 0.003 

tm tm 

Sb 

0.69 

— 

0*66 + 0*06 

— 

Ba 

675 

— 

705 + 48 

— 

a Using G- 

-2 as a standard. 




Best values given by Flanagan 

(17). 
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FIGURE CAPTIONS 


Figure 1 , 


Figure 2 . 

Figure 5 . 

Figure 4 . 
Figure 5 , 


Schematic diagram of target chamber and detection system. 

(A) vacuum fitting to allow target holder adjustment, (see 
enlargement inset), 

(B) target mounted on Aluminum holder, 

(C) quartz disc for focussing particle beam, 

(D) reference scale for sample positioning, 

(E) facuum fitting for Si(Li) detector connection to target 
chamber, o-rings indicated by solid circles, (see 
enlargement inset) , 

(F')beam stop. 

Characteristic X-ray spectra for USGS BCR-1 as a) 10$ graphit 
pellet and b) deposited on formvar film. 

Ionization cross sections for K and L X-ray production by 
4 MeV heavy charged particles. 

Gamma-yield for the 27 A1 (p,^) 2s Si reaction. 

Ratio of the Coulomb excitation cross sections for 4 MeV 


protons to alphas as a function of nucleus excitation level. 
Characteristic gamma-ray spectrum for USGS BCR-1. 



